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TRIGGERED STAR FORMATION AROUND MID-INFRARED BUBBLES IN G8.14+0.23 HII REGION 
L. K. Dewangan 1 , D. K. Ojha 1 , B. G. Anandarao 2 , S. K. Ghosh 3 , and S. Chakraborti 1 

ABSTRACT 

Mid-infrared (MIR) shells or bubbles around expanding Hll regions have received much attention 
due to their ability to initiate a new generation of star formation. We present multi-wavelength obser- 
vations around two bubbles associated with a southern massive star-forming (MSF) region G8. 14+0.23, 
to investigate the triggered star formation signature on the edges of the bubbles by the expansion of 
the Hll region. We have found observational signatures of the collected molecular and cold dust 
material along the bubbles and the 12 CO(J=3-2) velocity map reveals that the molecular gas in the 
bubbles is physically associated around the G8. 14+0.23 region. We have detected 244 young stellar 
objects (YSOs) in the region and about 37% of these YSOs occur in clusters. Interestingly, these YSO 
clusters are associated with the collected material on the edges of the bubbles. We have found good 
agreement between the dynamical age of the H II region and the kinematical time scale of bubbles 
(from the 12 CO(J=3-2) line data) with the fragmentation time of the accumulated molecular mate- 
rials to explain possible "collect-and-collapse" process around the G8. 14+0.23 region. However, one 
can not entirely rule out the possibility of triggered star formation by compression of the pre-existing 
dense clumps by the shock wave. We have also found two massive embedded YSOs (about 10 and 
22 M ) which are associated with the dense fragmented clump at the interface of the bubbles. We 
conclude that the expansion of the Hll region is also leading to the formation of these two young 
massive embedded YSOs in the G8. 14+0.23 region. 

Subject headings: dust, extinction - Hll regions - ISM: bubbles - ISM: individual objects (IRAS 
17599-2148) - stars: formation - stars: pre-main sequence 



1. INTRODUCTION 

Massive stars have a major impact on their surround- 
ing environment due to their energetic wind, UV ioniz- 
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ing radiation a nd an expanding Hll region (Zinnecker 
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In recent years, MIR shells or bub 
e expanding Hll regions have received 



much attention, because of th eir ability to initiate a new 
gener ation of star formati on (|Watson et al.| [2008, 2009, 
|2010[ ). IChurchwell et al.| fl2ofe[ |2007| ) identified aBout 
600 such shells or bubbles based on a circular morphology 
in Spitzer-GLIMPSE 8 fim emission. It is now believed 
that such regions are very promising to observationally 
stud y the conditions of seq uential /triggered star forma- 
tion (Thompson et al. |2012| and references therein). 

In this paper, we present a study of two MIR bub- 
bles, CN107 and CN109, from the catalogue of |Church-| 
well et al.| (|2007D around the G8. 14+0.23 region. The 



G8. 14+0.23 is an irregular Galactic ultra-compact (UC) 
H II region close to the IRAS 17599-2148 source , situ- 
ated at a distance of 4.2 kpc (|Kim fc Koo||200T|). The 



radio, sub-mm and molecular line observations dCodella 


et al. 


1994; Walsh et al.||1997l |1998l |Kim fc Koo 


2001, 


2003; 


Ojeda-May et al.|2002j [Thompson et al.|2006 


. The 



is also known as a bipolar blister type H H region due to a 
champagne flow ( [Kim fc KoopOOl] [2003] ) . The 6.7 GHz 
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Class II met hanol maser (Walsh et al. 1997 1998) and 
water maser (jCodella et al. 1994| ) were also detected close 



to the IRAS position. The 



site of massive star formation 



is often associated with methanol and water maser emis- 



sion. 



Kim fc Koo| ( |2001| ) found the velocity of ionized gas 
to be about 20.3 km s _i close to the IRAS position, using 
H76a recombination line profile. The velocity of molec- 
ular flow was o btained to be about 19.6 km s _1 using CS 
J= 2-1 lines by |Kim fc Koo (2003). These velocities are 
close to the values obtained by NH 3 (2,2) and (4,4) emis- 



sions to be about 19.3 km s (Churchwell et al. 



1990) 



and 20.3 km s 1 ( Cesaroni et al.|l992 ) , respectively. Kim 



fc Koo (2003) found irregular morphology and clumpy 
structure ol molecular cloud around G8. 14+0.23, using 
13 CO (J= 1-0) and CS (J= 2-1) line in t ensity maps. |Kim| 
fc Koo| ( |2003| and |Thompson et al.| fl2006|) also found 



a molecular and du st emission ridge around the region. 
|Okada et al. (2008) analyzed the MIR spectroscopic ob- 
servations using Spitzer-Infr&red Spectrograph (IRS) for 
14 Galactic star- forming regions, including G8. 14+0.23, 
in the wavelength range of 20 - 36.5 fim. They argued the 
presence of the shocked and ionized gas associated with 
the G8. 14+0.23 region, based on the detection of [Sill], 
[Fe n] and [S ill] emission lines and derived the abun- 
dance ratio of two ions against the solar abundance (as) 
as 0.060, 0.158 and 0.020 for (Fe+/Si+) as , (Si+/S 2+ ) as 
and (Fe 2+ /S 2+ ) as respectively. 

In recent years, multi- wavelength observations are uti- 
lized to trace out different components associated with 
the star-forming regions viz. ionized, shocked, cold dust 
and molecular £ as, as well as embedded young stellar 
populations. In this paper, we present multi- wavelength 
observations to investigate the triggered star formation 
around the MIR bubbles associated with a southern MSF 
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region G8. 14+0.23. This region is not well explored in 
infrared wavelengths, especially the identification of em- 
bedded populations around the region. Our study also 
includes the search for the infrared counterpart (s) of em- 
bedded massive YSOs in their earlier phases and to un- 
derstand their formation processes. 

In Section [2j we introduce the archival data and data 
reduction procedures used for the present study. In Sec- 
tion[3j we examine the morphology of the G8. 14+0.23 re- 
gion in different wavelengths and the interaction of mas- 
sive stars with its environment using various Spitzer MIR 
ratio maps. We also describe the selection of young pop- 
ulation, their distribution around the G8. 14+0.23 region 
and discuss the triggered star formation scenario on the 
edge of the bubbles in this section. In Section [4j we 
summarize our conclusions. 



2. AVAILABLE DATA AROUND G8. 14+0.23 AND DATA 
REDUCTION 

Archival deep near-infrared (NIR) JHK S images and 
catalogue of G8. 14+0.23 region were obtained from the 
UKIDSS 6 th archival data release (U KIDSSDR6plus) 
of th e Galactic Plane Survey (GPS) ( jLawrence et al 
'2007b. UKIDSS observations were made using the 



UK1RT Wide Field Camera (WFCAM; ICasali et al. 



2007) and fluxes were calibrated using Tw o Micron All 
S^TBurvey (2MASS; [Skrutskie et aLp OQ6). The details 
of basic dat a reduction and cali bration procedures ar e 
described in |Dye et al.] ( |2006|) and |Hodgkin et al .M2009| ), 
respectively. Magnitudes ol bright stars (J ^ ll.o mag, 
H $C 11.5 mag and K s ^ 10.5 mag) were obtained from 
the 2MASS, due to saturation of UKIDSS bright sources. 
Only those sources are selected for the study which have 
photometric magnitude error of 0.1 and less in each band. 

We obtained narrow band molecular hydrogen (H2; 
2.12 /im; 1 - S(lV) ima ging data from UWISH2 sur- 
vey ( Froebrich et al.|2011[). We follow ed a similar proce- 
dure ai~cIe^clTEe^n5y^^mca^ ( |2011[ ) to obtain the final 
continuum-subtracted H2 image using GPS K s image. 

The Spitzer Space Telescope Infrared Array Camera 
(IRAC (Chi ( 3.6 /im), Ch2 (4.5 am), Ch3 (5.8 /im) and 



Ch4 (8.0 /im); |Fazio et aLl|2004[ )) and Multiband Imag 
ing Photometer (MIPS (2417m); |Rieke et al] fl2004| )) 
archival images were obtained around the G8.14+0.z3 
region from the "Galactic Legacy Infrared Mid-Plane 
Survey Extraordinaire" (GLIMPSE) and "A 24 and 70 
Micron Survey of the Inner Galactic Disk with MIPS" 
(MIPSGAL) surveys. MIPSGAL 24 /im image is satu- 
rated close to the IRAS 17599-2148 position. We per- 
formed aperture photometry on all the GLIMPSE im- 
ages (plate scale of 0.6 arcsec/pixel) using a 2.4 arcsec 
aperture and a sky annulus from 2.4 to 7.3 arcsec using 
IRAF 1 . The photometry is calibrated using zero mag- 
nitudes including aperture corrections, 18.5931 (Chi), 
18.0895 (Ch2), 17.4899 (Ch3) and 16.6997 (Ch4), ob- 
tained from IRAC Instrument Handbook (Version 1.0, 
February 2010). 

We obtained 20 cm radio continuum map (resolution 
~ 6 arcsec) from Very Large Array (VLA) Multi- Array 
Galactic Plane Imaging Survey (MAGPIS) jHelfand et] 
|al.|2006[ ) to trace the ionized region around G8. 14+0.23. 

1 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, USA 



The molecular 12 CO(J=3-2) (rest frequency 345.7959899 
GHz) spectral line public processed archival data was 
also utilized in the present work. The CO observations 
were taken on 11 June 2009 at the 15 m James Clerk 
Maxwell Telescope (JCMT) using the HARP array. The 
JCMT archival reduced fits cube was used to extract the 
velocity information around the region using the Astro- 
nomical Image Processing System 2 (AIPS) package. The 
fits cube was loaded into AIPS and rotated to make the 
frequency axis the first one. An alternate velocity axis 
was defined from the frequency using the rest frame fre- 
quency of the CO J=3-2 transition. The data were Han- 
ning smoothed over 476 m s _1 to decrease the noise per 
channel. It was then integrated over the velocity range 
10 to 30 km s _1 to obtain the zeroth moment (column 
density), first moment (mean velocity) and second mo- 
ment (velocity dispersion). Archival BOLO CAM 1.1 mm 
jAguirre et"adT|2Qll| ) and SCUBA 850 urn ( |Di Francesco 
et al. 2D08| ) images were also used in the present work. 
The effective FWHM Gaussian beam sizes of final 1.1 
mm and 850 /am images were 33 arcsec and 22.9 arcsec, 
respectively. 



3. RESULTS AND DISCUSSION 



3.1. 



Morphology of the G8. 14+0.23 region in different 
wavelengths 

Figure represents the RGB color composite image 
using GLIMPSE (8.0 /im (red) & 4.5 jam (green)) and 
UKIDSS K s (blue) of a region (~ 6.6 x 5.9 arcmin 2 ) 
around G8. 14+0.23. The 8 /im band contains the two 
strongest polycyclic aromatic hydrocarbon (PAH) fea- 
tures at 7.7 /im and 8.6 /im, which are excited in the 
photodissociation region (or photon-dominated region, 
or PDR). The PDRs are the interface between neu- 
tral & molecular hydrogen and traced by PAH emis- 
sions. Figure [lb shows the 3 color composite image using 
UKIDSS J (blue), H (green) and K s (red) band in log 
scale. The positions of IRAS 17599-2148 (*), UC H11 
region (O), methanol maser (+) and water maser (x) 
are marked in the figure. The coordinates of UC H II re- 
gion, methanol maser and water maser were taken from 



Kimfc Koo (20011), Walsh et al. (1998) and Codella et al. 



(|1994[ ), respectively. Figure[l|i displays the two extended 
1V11K bubbles CN107 and CN109 prominently around the 
G8. 14+0.23 region. It is to be noted that these bubble 
structures are not visible in any of the UKIDSS NIR 
images (see Figure jib), but dark regions are seen corre- 
sponding to the bubble structures in Figure [TJd. The two 
infrared sources (designated as IRS 1 and IRiS2) are seen 
close to the IRAS pos ition a nd are also marked in Fig- 
ure [l^i (see subsection 3.3.1 for more details). Figure [2] 
shows color composite image made using MIPSGAL 2T 
jam (red), GLIMPSE 8 /im (green), and 3.6 jam (blue) 
images, overlaid by BOLOCAM 1.1 mm and SCUBA 850 
/am emission by dotted white and solid yellow contours, 
respectively. MIPS 24 /am image is saturated near to the 
IRAS position. MAGPIS 20 cm radio continuum emis- 
sion is also overlaid in Figure [2] by black contours. The 
cold dust emission at 850 jam is very dense and prominent 
at the interface of the two bubbles. It is to be noted that 
the peaks of cold dust and ionized gas are concentrated 

2 http://www.aips.nrao.edu/ 
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Fig. 1. 



a) 3 



6.6 



5.9 arcmin 



color composite image (size 

central coordinates: a 2 ooo = 18 h 02™ 56 s . 3, 5 2 ooo = -21° 48 12 .6) 
of the G8. 14+0.23 region, using Spitzer-GLIMPSE images at 8.0 /im 
(red), 4.5 fim (green) and UKIDSS K s (blue) in log scale. The positions 
of the two sources (IRS 1 & IRS 2) are marked by blue circles and 
labeled on the image. The scale bar on the top left shows a size of 1 pc 
at the distance of 4.2 kpc. The white dash ed box is shown as zo omed- 
in view in Fig. [8] Two bubbles from the [Churchwell et al. ( 2007} are 
also labeled on the image as CN 107 and CN 109. b) Color composite 
image using UKIDSS J (blue), H (green) and K s (red) in log scale. The 
positions of IRAS 17599-2148 (*), UC Hn region (O), methanol maser 
(+) and water maser (x) are marked in the figure. 

near the IRAS position and interface of the bubbles. The 
24 /im and 20 cm images trace the warm dust and ion- 
ized gas in the region, respectively. It is obvious from 
r lgure [2] that the PDR region (traced by 8 /im) encloses 
24 /xm and 20 cm radio emissions inside the bubbles and 
indicates the presence of dust in and arou nd the H II re- 
gion (see for example | Wat son et al-|[2008| for N10, N21, 
and N49 bubbles). 

3.2. IRAC Ratio Maps and the Collected material 

In recent years, Spitzer-IRAC bands and ratio maps are 
utilized to study the interacti on of massive s tars with its 
immediate environment ( Povich et al.||2007 ). The IRAC 
bands contain a number of prominent atomic and molec- 
ular lines /fe atures such as H2 line s in all channels (see 
Table 1 from |Smith fcRos en 2005), Bra 4.05 /im (Ch2), 
Fell 5.34 /im (Ch3), Aril 6.9 9 urn and Arm 8.99 /im 
(Ch4) (see Reach et al. 2006). It is to be noted that 
the Spitzer-IRAC bands, Chi, Ch3 and Ch4, contain the 
PAH features at 3.3, 6.2, 7.7 and 8.6 /im, whereas Ch2 
(4.5 /im) does not include any PAH features. Therefore, 
IRAC ratio (Ch4/Ch2, Ch3/Ch2 and Chl/Ch2) maps 
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Fig. 2. — Figure represents 20 cm contours in black color around 
the G8. 14+0.23 region obtained from the MAGPIS survey, overlaid 
on a color composite image made using the 24 /im (red), 8 /im 
(green), and 3.6 /im (blue) images. The 20 cm contour levels are 
10, 25, 40, 55, 70, 85 and 95 % of the peak value i.e. 0.23 Jy/ 
beam. BOLOCAM 1.1 mm and SCUBA 850 /im emissions are 
also shown by dotted white and solid yellow contours with 10, 25, 
40, 55, 70, 85 and 95 % of the peak valu e i.e. 3.55 and 8.23 Jy/ 
beam, respectively. Two bubbles from the Church well et al.] ( 2007 ) 
are also labeled on the image as CN 107 and (JIN 1UU. The marked 
symbols are similar to as shown in Fig. fl] (see text for more details). 
MIPS 24 /im image is saturated near to the IRAS position. 

are being used to trace out the PAH features in MSF re- 
gions (e.g. | Povich et al.|2007||Watson et al.|2QQ8| [Kumar 
& Anandarao 20lU| due to UV radiation from massive 
star(s). In order to make the ratio maps, we generate 
residual frames for each band removing point sources by 
choosing an extended aperture (1 2.2 arcsec) and a l arger 
sky annulus (14.6 - 24.4 ar csec; pleach et al. 2005) in 
IRAF/DAOPHOT software flStetson|1987| . These resid- 
ual frames are then subjected to median filtering with a 
width of 10 pixels and s moothing by 3 x 3 pixels using 
the boxcar algorithm (e.g. Povich et al.|2007 ). Figures |3p 
and [3b represent the IRAC ratio maps, Ch3/Ch2 and 
Ch4/Ch2 around the G8. 14+0.23 region, respectively. 
The ratio contours are also overlaid on the ratio maps 
for better clarity and insight (see Fig. |3j) . Both the ra- 
tio maps clearly trace the prominent PAH emissions and 
subsequently extent of PDRs in the region. 

Figure p] shows the IRAC Ch2/Chl ratio map over- 
laid with (Jh3/Ch2 ratio and MAGPIS 20 cm radio con- 
tours. The bright region traced by the Ch2/Chl ratio 
map is coincident with the 20 cm radio continuum emis- 
sion, possibly due to presence of the Bra (4.05 /im) fea- 
ture in Ch2 band around the Hll region. One can also 
notice very faint diffuse emission along the edges of the 
bubbles away from the peak of 20 cm radio emission, 
possib ly due to the molecular hydrogen feature in Ch2 



band. |Ojeda-May et al.| (2002) detected two sources in 
the 6 cm radio observations around the G8. 14+0.23 re- 
gion. The positions of these sources are also shown in 
Figure [4] Figures [5^i, [5J3 & [5b show the grey-scale im- 
ages in SCUBA 850 /im, BOLOCAM 1.1 mm and JCMT 
CO 3-2 (345.79599 GHz) around the G8.14+0.23 region 
respectively. Figure [5] exhibits the evidence of collected 
material along the bubbles traced by dust continuum 
(SCUBA 850 /im & BOLOCAM 1.1 mm) and molec- 
ular gas (JCMT CO 3-2 at 345.79599 GHz) emission, 
which is also indicated by red arrows. These arrows are 
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FIG. 3. — a) Figure exhibits the GLIMPSE Ch3/Ch2 ratio map of 
the G8. 14+0. 23 region (similar area as shown in Fig. fTl. The ratio 
Ch3/Ch2 value is found to be "7—8" and "8.5—9" in the bubble interior 
and for the brightest part of the PDR respectively. The Ch3/Ch2 
ratio contours are also overlaid on the image with a level of 7.7, a 
representative value between "7—8". b) The Ch4/Ch2 ratio map of the 
region is shown here. The ratio Ch4/Ch2 value for the brightest part 
of the PDR and the interior is "30 34" and "23 25" respectively. The 
Ch4/Ch2 ratio contours are also overlaid on the image with a level of 
24.8, a representative value between "23 25". The marked symbols are 
similar to as shown in Fig. ^ 
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Fig. 4. — Figure exhibits the Ch2/Chl ratio image of the 
G8. 14+0.23 region. The Ch3/Ch2 ratio contours in violet are also 
overlaid on the image with a level of 7.7. Contours in red color 
represent the 20 cm radio continuum emission with similar levels 
as shown in Fig. [2] around th e region. Square symbol s represent 
the 6 cm radio detections by |Ojeda-May et al. ( 2002 ) and other 
marked symbols are similar to as shown m Fig! ^ 
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Fig. 5. The grey-scale image in SCUBA 850 fim (a), BOLO- 
CAM 1.1 mm (b) and JCMT CO 3-2 (c) around the G8. 14+0.23 
region, respectively. Red arrows show the detected faint features 
around the G8. 14+0.23 region and other marked symbols are sim- 
ilar to as shown in Fig. [I] 



marked to show the associated faint features above the 
4-sigma noise level for each image. We further utilized 
JCMT public processed CO 3-2 cube data to observa- 
tionally check the association of the molecular material 
with the expanding H II region seen as MIR bubbles. We 
calculated the mean H2 number density near the H II re- 
gion using 12 CO zeroth moment map. We derived the 
column density using the formula 7Vh 2 (cm -2 ) = X x 
Wcq, w here X = 6 x 10 20 cm -2 K -1 km -1 s (see Ji et 
al.|2Q12[ ) and W C o = 48.08 K km s -1 from our map. The 



mean H2 number density is obtained to be 3575.7 cm 
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Fig. 6. — Figure shows the 12 CO(J=3-2) channel maps in 1 km 
s — 1 steps with velocity interval between 17 to 25 km s _1 around 
the G8. 14+0.23 region. The grey color scale bar at the top is 
presented in unit of K in the range of 4 to 32 K and the contours 
are plotted for 6.4, 12.8 and 25.7 K, from outer to inner side. The 
bubbles CN 107 and CN 109 are also labeled in the top left panel. 



using the relation 7Vh 2 (cm 2 )/L (cm), where L is the 
molecular core size of about 8.07 x 10 18 cm (~ 2.6 pc) 
near the H II region. Figure [6] exhibits molecular emission 
around the region between 17 and 25 km s _1 velocity in- 
terval in steps of 1 km s _1 . The bubbles (CN 107 and 
CN 109) are clearly traced by molecular emission in the 
velocity range of 19 to 21 km s _1 (see Fig. |6|. We also 
estimated velocity dispersion of about 4 km s _1 (using 
second moment CO map) around the region, assuming 
a smooth rotation curve for the Galaxy with typical gas 
dispersion veloc ity of about 10 km s _1 ( |Stark fc Brand| 
|1989| )Silk||1997| ). The consideration of velocity ranges 
ol molecular gas in which bubbles are traced (see Fig. [6]) 
and the velocity dispersion value, imply that the molecu- 
lar gas in these bubbles are physically associated. These 
velocity ranges of molecular gas are also compatible with 
the ionized gas velocity (~ 20.3 km s _1 ) obtained by the 
H76a r ecombination line profile study around the Hll 
region ( Kim fc Koop OOl), which shows the physical as- 
sociation of the molecular material and the Hll region. 
The evidence of collected material along the bubbles is 
further confirmed by the detection of the H2 emission (see 
Figure [7]). Figure [^represents the continuum-subtracted 
H2 image at 2.12 jam and reveals that the H2 emission 
surrounds the Hll region along bubbles, forming a PDR 
region, which may be collected due to the shock. In brief, 
the PAH emission, cold dust emission, molecular CO gas 
and shocked H2 emissions are coincident along the bub- 
bles. 

Figure [8] exhibits the zoomed-in color composite image 
using IRAC/GLIMPSE images (red: 8.0 /im, green: 5.8 
jam and blue: 4.5 jam) around the G8. 14+0.23 region. 
The slit positions of Spitzer-IRS spectrograph at 8 ob- 
served locations around the G8. 14+0.23 region are also 
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Right Ascension (J2000) 

Fig. 7. — Figure represents the inverted grey-scale image of the 
continuum-subtracted H2 image at 2.12 /im. White arrows indicate 
the detected H2 emission along the bubbles (as seen in GLIMPSE 
images) around the G8. 14+0.23 region. The other marked symbols 
are similar to as shown in Fig. ^ The continuum-subtracted H2 
image is processed to median filtering with a width of 4 pixels and 
smoothened by 4 x 4 pixels using boxcar algorithm to trace out 
the faint features in the image. 
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Fig. 8. — Zoomed-in 3 color composite image using IRAC 4.5 /im 
(blue), 5.8 /im (green) and 8.0 /im (red) around the G8. 14+0.23 
region (size ~ 2.36 x 1.26 arcmin 2 ; as shown by a dashed box 
in Fig. flj. Star symbols represent 8 positions of Spitzer-IRS slit 
around The G8. 14+0.23 region and also labeled as to 7. SCUBA 
850 /im emission is also shown by dotted yellow contours with 
similar levels as shown in Fig. [2] The positions of IRS 1 & IRS 2 
sources are marked by red circles and labeled in the image. 



marked in the Figure [8] with star symbols (see |Okada et 
|aL ||2008| ). |Okada et aL| (|2'0Q8) also listed the line inten- 
sities ol detected emission lines ([Arm], 22 /im; [Fein], 
23 /im; [Fell], 26 /im; H 2 0-0 S(0), 28 /im; [S ill], 33 
/im; [Sill], 35 /im and [Nelll], 36 jam) at the 8 observed 
positions (similar labels marked as to 7 in Figure [8 
around the G8. 14+0.23 region. The presence of bot 
PAH and fine structure line emissions around the region 
clearly suggest, the presence of an ionizing source with 
UV radiation close to radio and dust continuum peaks. 

3.3. Photometric analysis of point-like sources around 
the G8. 14+0. 23 region 

In order to trace ongoing star formation activity in the 
G8. 14+0.23 region, we have identified YSOs using NIR 
and GLIMPSE data. 

3.3.1. Identification of YSOs 



We have used Gutermuth et al. (2009) criteria based 
on four IRAC bands to identity 1 SOs and various 
possible contaminants (e.g. broad-line AGNs, PAH- 
emitting galaxies, shocked emission blobs/knots and 
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PAH-emission-contaminated apertures). These YSOs 
are further classified into different evolutionary stages 
(i.e. Class I, Class II, Class III and photospheres) using 
slopes of the IRAC spectral energy distribution (SED). 
Fig. [9^ shows the IRAC color-color ([3.6]-[4.5] vs [5.8]- 
[8.0] jdiagram for all the identified sources. We find 14 
YSOs (6 Class 0/1; 8 Class II), 2 Class III, 123 photo- 
spheres and 8 contaminants in the G8. 14+0.23 region. 
Two sources close to the peak of the dense molecular gas 
and the dust clump at the interface of the bubbles, IRS 1 
and IRS 2, are also labeled in Fig. [9^i. The position of 
IRS 2 is shown by a star symbol using 3.6 /im mag as an 
upper limit (see Table pi). It is to be noted that IRS 2 is 
close (about 2 arcsec) to the detected methanol maser. 
The details of the YSO classi fication can be found in 
|Dewangan fc Anandarao| ( |2011[ ) . We have also applied 
criteria to i dentity Y SUs as suggested by |Hartmann et aL 
(I2005D and IGetman et all (120071) for those sources that 
are detected in three IRAC/ GLIMPSE bands, but not 
in the 8.0 /im band. We identify 12 more YSOs (2 Class 
0/1; 10 Class II) through color-color diagram using three 
GLIMPSE bands in the region (see Fig. ^p). 

UKIDSS NIR JHK S photometry is used to identify 
more YSOs in the region. NIR color-color diagram (CC- 
D: H-K s /J-H) is shown in Fig.[9fc for all the sources de- 
tected in J, H and K s bands. We have divided the CC 
diagr am into three regions namely "F", "T", an d "P" 
(e.g., |Sugitani et"aI1|2002| |Ojha et al.)|2004a|b[ ). The 
"F" sources are located between the reddening bands of 
the main-sequence and giant stars. The sources located 
within the "T" region along the T Tauri locus are T 
Tauri-like sources (Class II objects) with large NIR ex- 
cess. "P" sources are those located in the region redward 
of region "T", and are most likely Class I objects (pro- 
tostellar objects). We have identified 159 Class II and 
30 Class I sources located in the "T" and "P" regions 
in the Fig. [9^, respectively. The identified Class I and 
Class II sources are shown by red circles and blue trian- 
gles in Fig. |9j3. We have estimated the visual e xtinction 
Lsing the CC diagram and following the |Indebe-| 
extinction law, for those sources that 

The 



touw et al 
tall within 



using tne 00 
et al. (120051) e 
Ithrn tEPT" 



2010) 



region (see Ojha et al 
average value comes out to be Ay ^8.2 mag, which is 
in good agreement with that estimated (Ay ~ 8.8 mag) 
by |Rowles fc Froebrich| Q2009| ) using 2MASS data for the 
G8. 14+0.23 region. 

GLIMPSE 3.6 and 4.5 /im bands are more sensitive 
for point sources than GLIMPSE 5.8 and 8.0 fim images. 
Therefore, a larger number of YSOs can be identified us- 
ing a combination of NIR (JHK S ) with GLIMPSE 3.6 
and 4.5 /im (i.e. NIR-IRAC) photometry, where sources 
are not detected in IRAC 5.8 and/or 8.0 /im band ( |Guter-| 
muth et al.||2009 |) . We followed the criteria given by 



Gutermuth et 8lL\ ( |2009[ ) to identify YSOs using H, K s 
3.6 and 4.5 /im data. We have found 29 more YSOs 
(26 Class II and 3 Class I) using NIR-IRAC data (see 
Fig.gl). 

Finally, we have obtained a total of 244 YSOs (203 
Class II and 41 Class I) using NIR and GLIMPSE data 
in the region. 

3.3.2. YSOs surface density map 

To study the spatial distribution of YSOs, we gener- 
ated the surface density map of all YSOs using a 5 arcsec 




.;.;/ 









FlG. 9. — a) Color-color diagram (CC-D) using the Spitzer-IRAC 
four bands for all the sources identified within the region shown in 
Fig. ^ The extinction vector for Ak = 5 mag is shown by an arrow, 
using average extinction law from Flaherty et al. (2007). The dots in 
grey color around the centre (0,0; locate the stars with only photo- 
spheric emissions. The open squares (black), open triangles (blue) and 
open circles (red) represent Class III, Class II and Class 0/1 sources 
respectively, classified using the olirac criteria. The "x" symbols in 
magenta color show the identified PAH-emission-contaminated aper- 
tures in the region. The source IRS 2 is marked by a red star symbol 
and is also labeled along with the source IRS 1 in the diagram (see 
the text), b) CC-D of the sources detected in three IRAC/GLIMPSE 
bands, except 8.0 /im. The open squares (black), open triangles (blue) 
and open circles (red) represent Class III, Class II and Class 0/1 sources 
respective ly, cl assified following the criteria suggested by |Hartmann et] 
[aI1 ( [2005t and jGetman et al.|p007). c) UKIDSS JHK S CC-D for the 
G8.14+0.23 region. UKIDSS JN1R (JHK S ) magnitudes are calibrated 
in 2MASS system. The unreddened dwarf and giants loci are shown 
by red and magenta curves ( jBessell &; Brett]|l988| ) respectively. The 
long-dashed straight lines across the CC-D are the extinction vectors 
(Ak = 3 mag) using Indebetouw et al. (2005J extinction law. The ex- 
tinction vector Ak = 1 mag is also shown m the diagram. Classical T 
Tauri (CTTS) locus (in California Institute of Technology (CIT) sys- 
tem) jMeyer et al.|l997| ) is shown by a dashed-dotted line. The loci of 
unreddened dwarf (Bessell & Brett (BB) system), giant (BB-system) 
and CTTS (CIT system) are c onverted int o 2MASS system using trans- 
formation equations given by Carpenter (2001). Red circles and open 
blue triangles represent the Class I and class II YSOs. The CC-D is 
classified into three regions namely "F" , "T" , and "P" (see text for 
details), d) Figure shows the de-reddened [K s - [3.6]]o vs [[3.6] - [4.5]]o 
CC-D using NIR and GLIMPSE data. The selected region shown by 
the solid lines represents the location of YSOs. The extinction vector 
for Ak = 2 mag is shown by an arrow, calculated using the average 
extinction law from Flaherty et al. (2007). Open red circles and open 
represent Class I and UTi 



blue triangles represent < 



Class II sources respectively. 



grid size, fol l owing the same procedure as given in |Guter- 
muth et al. (2009). The surface density map of YSOs is 
constructed using 6 nearest-neighbour (NN) YSOs for 



each grid point. Fig. 
of all identified YSOs 



i shows the spatial distribution 
ass I and Class II) in the region. 
The contours of YSO surface density and Ch3/Ch2 ra- 
tio map are also overlaid on the diagram. The levels 
of YSO surface density contours are 6 (2.5a), 8 (3.3a), 
11 (4.6cr) and 16 (6.7a) YSOs/pc 2 , increasing from the 
outer to the inner region. We have also calculated the 
empirical cumulative distribution (ECD) as a function of 
NN distance to identify the clustered YSOs in the region. 
Using the ECD, we estimate the distance of inflection d c 
= 0.617 pc (0.0084 degrees at 4.2 kpc) f or the region 
for a surface den sity of 6 YSOs/pc 2 (see |Dewangan fc| 
Anandarao 2011[ for details of d c and ECD). We 
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TABLE 1 

NIR and Spitzer IRAC/GLIMPSE photometric magnitudes are listed for selected clustered YSOs having less than or equal to 

d c VALUE CLOSE TO THE SUB-MM PEAK (SEE THE TEXT) . SPECTRAL INDEX (a IR Ac) WAS OBTAINED BY FITTING OF IRAC WAVELENGTHS 

(sources having atleast 3 IRAC wavelengths magnitude). 



Source 


RA[2000] 


DEC[2000] 


J (mag) 


H (mag) 


K s (mag) 


Chi (mag) 


Ch2 (mag) 


Ch3 (mag) 


Ch4 (mag) 


OLIRAC 


IRS 1 


18 


03:00.94 


-21 


48:07.3 






13.44±0.01 


8.79±0.01 


7.29±0.01 


6.07±0.03 


4.67±0.04 


1.85 


IRS 2 


18 


03:00.59 


-21 


48:10.5 








< 11.78 


10.39±0.10 


7.68±0.05 


6.13±0.07 


3.96 


IRS 3 


18 


03:00.36 


-21 


48:23.5 


15.97±0.01 


13.76±0.00 


12.50±0.00 


11.44±0.08 


10.83±0.09 








IRS 4 


18 


03:01.12 


-21 


48:01.9 


16.63±0.02 


14.98±0.01 


13.54±0.01 


11.10±0.14 


9.43±0.10 








IRS 5 


18 


03:02.90 


-21 


48:23.1 






15.27±0.03 


11.54±0.02 


9.92±0.02 


8.71±0.04 


7.44±0.06 


1.81 


IRS 6 


18 


03:03.70 


-21 


48:45.8 


17.94±0.05 


15.87±0.02 


14.51±0.01 


13.14±0.05 


12.59±0.07 








IRS 7 


18 


03:04.33 


-21 


47:46.1 


15.20±0.01 


14.39±0.01 


13.88±0.01 


13.44±0.08 


13.05±0.10 








IRS 8 


18 


03:04.39 


-21 


48:31.1 






15.72±0.04 


11.27±0.01 


9.47±0.01 


8.55±0.01 


7.79±0.02 


1.03 



that 37% (91 out of 244 YSOs) YSOs are present in clus- 
ters. Figure [TO^ i reveals that the distribution of YSOs 
is mostly concentrated in the North and East regions, 
having peak density of about 20 YSOs/pc 2 , while YSO 
density of about 8 (3a) YSOs/pc 2 is also seen around 
the West and South-East regions (see Figure [Tok ). These 
values of YSO surface density is much lowertnan usual 
value s (about 30-60 YSOs/pc 2 ) found in recent works 
(see |Gutermuth et al |[2009| |2011[ ) in MSF regions. Also 
the percentage of YSOs in clusters is low (only 37%) as 
compared to recent reports (about 60-80%) in MSF re- 



gions (see IGutermuth et al.||2009| |2011 



Chavarria et al. 

|2008[ ). It possibly supports the youth of the region and 
a site of further ongoing star formation. It is also to be 
noted that the YSO surface density is associated with 
the PDR region, on the edge of bubbles. The correlation 
of cold dust, molecular gas, ionized gas and YSO surface 
density is shown in Figure [To|d. The association of YSOs 
with the collected materials around the region further re- 
veals the ongoing star formation on the edge of bubbles. 
In addition, we have checked the possibility of intrinsi- 
cally "red sources" contami nation, such as AGB sta rs in 
our YSO sample. Recently, Robitaille et al. (2008) pre- 
pared an extensive catalogue of such red sources based 
on the Spitzer GLIMPSE and MIPSGAL surveys. They 
showed that two classes of sources are well separated in 
the [8.0 - 24.0] color space such that Y SOs are redder 
than AGB stars in this space (see also | Whitney et al.| 
2008|). Since, MIPS 24 /im data is saturated lor the 



G8. 14+0.23 region, we have therefore used criteria based 
on the IRAC magnitude (4.5 /im) and color space ([4.5 
- 8.0]), to estimate the AGB contamination. It is be- 
lieved tha t the presence of AGB stars in th e cluster is un- 
likely (see |Dewangan fc An andarao 20lT| and references 
therein). Therefore, we have applied red source criteria 
for only those YSOs which are cluster members with less 
than or equal to d c value (see Fig. 10 1). We find that 
GLIMPSE YSOs may be contaminate! by AGB stars u p 
to about 36% (5 out of 14 YSOs; see subsection 3.3.1). 



3.3.3. SED modeling of clustered YSOs 

In this subsection, we present SED modeling of some 
selected YSOs to derive their various ph ysical param- 
eters using an on-line SED modeling tool (Robitaille et| 



al.|2006[ [2007 ) . Clustering analysis of identified YSOs re- 
veals that there is a grouping of YSOs (~ 20 YSOs/pc 2 ) 
close to the peak of the dense molecular gas and the dust 
clump in the region (see Fig. [lofr ). We have selected 8 
YSOs (designated as IRS 1 IRS 8) from this cluster 




18 h 05 m 1 s 5 s s 02 m 55 s 50 s 45 s 

) 1 PC 



. [4.5 pj 




s 02 m 55 s 50 s 

Right Ascension (J2000) 

FlG. 10. — a) Figure represents the spatial distribution of all identi- 
fied YSOs in the G8. 14+0.23 region. The YSO surface density contours 
are plotted for 6, 8, 11 and 16 YSOs/pc 2 , from outer to inner side (see 
text for details). The open circles and open triangles show the Class I 
and Class II sources respectively. The YSOs identified using four IRAC, 
three IRAC, NIR-IRAC and NIR data are shown by red, orange, green 
and blue colors respectively. The Ch3/Ch2 ratio contours in violet 
color are also overlaid on the image with a similar level as in Fig. ^] 
and other marked symbols are similar to as shown in Fig. ^ b) The 
schematic represents the spatial distribution of YSOs, dust emission, 
molecular gas and ionized gas in the G8. 14+0.23 region. GLIMPSE 
4.5 /im image is overl aid w ith contours of YSO surface density (blue; 
similar levels as in Fig.flO^), SCUBA 850 /im (yellow), MAGPIS 20 cm 
(cyan) and JCMT CO 3^2 (violet red) data. The CO 3-2 contour levels 
are 20, 29, 40, 55, 70 and 80 % of the peak i.e. 139.551 K km s _1 . The 
contour levels of SCUBA 850 /im and MAGPIS 20 cm are similar to 
as shown in Fig. [2] The selected clustered YSOs (see subsection 3.3.3) 
are also marked (red circles) and labeled as 1,....,8 on the image (see 
Table [TJ. 
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close to the dense clump, which are detected at least upto 
GLIMPSE 4.5 /im or longer wavelengths for SED mod- 
eling. NIR and Spitzer IRAC/GLIMPSE photometric 
magnitudes for these selected YSOs are listed in Table [I] 
along with IRAC spectral indices (cvirac) an d are also 
labeled in Fig. [To}}. IRAC spectral indices were calcu- 
lated using a least squares fit to the IRAC flux points 
in a log(A) versus log(AFA) diagram for thos e sources 



that are detected in at least 3 IRAC bands (see Dewangan 
'fc Anandarao||2QTTj for details). We have not included 



the MIPS 24 /im magnitude for SED modeling because 
the image is saturated close to the IRAS position. The 
SED model tool requires a minimum of three data points 
with good quality as well as the distance to the source 
and visual extinction value. A distance range of 3.7 to 
4.7 kpc and the visual extinction in the range from 8 to 
25 mag are used as input parameters in SED modeling 
tool for each source to constraint the SEDs. These mod- 
els assume an accretion scenario with a central source 
associated with rotationally flattened infalling envelope, 
bipolar cavities, and a flared accretion disk, all under ra- 
diative equilibrium. The model grid consists of 20,000 
models of two-dimensional Monte Carlo simulations of 
radiation transfer with 10 inclination angles, resulting in 
a total of 200,000 SED models. The grid of SED models 
covers the mass range from 0.1 to 50 M . Only those 
models are selected that satisfy the criterion \ 2 ~ Xbest 
< 3, where x 2 is taken per data point. The plots of SED 
fitted models are shown in Fig. [TT] for all the selected 
sources. The weighted mean values of the physical pa- 
rameters (age, mass, temperature, luminosity, envelope 
accretion rate (M env ) and disk accretion rate (Mdisk)) 
along with the standard deviations derived from the SED 
modeling for all the selected sources are given in Table [2] 
The table also contains the model-derived weighted mean 
values of Ay with standard deviations and the degener- 
acy of the models (i.e., the number of models that satisfy 
the x 2 criterion as mentioned above). The derived SED 
model parameters show that the average value of mass, 
age and Ay of these clustered YSOs are about 7.8 M , 
0.12 Myr and 15.5 mag, respectively. It is interesting 
to note that the embedded IRS 1 and IRS 2 sources are 
young and massive protostars. Finally, the SED mod- 
eling results favour the ongoing star formation around 
the region with detection of young YSOs as well as some 
massive candidates in their early phase of formation. 

3.4. Star formation scenario 

We have found evidence of collected material along 
the bubbles and also ongoing formation of YSOs on the 
edge of the bubbles. The morphology and distribution 
of YSOs suggest that the G8. 14+0.23 region is a site of 
star formation possibly triggered by the expansion of the 
Hll region. In recent years, the triggered star formation 
process, especially "collect and collapse" mechanism has 
been studied extensively on the edges of many Hll re- 
gions such as Sh 2-104, RCW 79, Sh 2-212, RCW 120, 



Sh 2-217 ( 


Deharveng et al.|2003, 


2008 


2009 Zavagnoet 


al.|2006 2 


010 Brand et al.| 2011 


). In order to check the 



around the G8. 14+0.23 region, we have calculated the 
dynamical age (td yn ) of the H 1 1 region and compar ed it 
with an analytical model by | Whitworth et ah] ([1994 ) . We 




Fig. 11. — The SED plots of all the selected sources including 
IRS 1 & IRS 2 are shown here. Filled circles are observed fluxes 
of good quality (with filled triangles as upper limits) taken from 
the archives or published literature (see text for more details) and 
the curves show the fitted model with criteria % 2 - Xbest ^ 3- ^ ne 
thin black curve corresponds to the best fitting model. The dashed 
curves represent photospheric contributions. 

have estimated the age of the Hll re gion at a given ra- 



dius R, using the following equation (|Dyson & Williams 
1980] ): 

4#, 



^dyn — 



7c a 



7/4 



(i) 



where c s is the isothermal sound velocity in the ion- 
ized gas (c s = 10 km s _1 ) and R s is the radius of the 
Stromgren sphere, given by R s = (3 N^^Tra^ae) 1 / 3 , 
where the radiative recombina tion coefficie nt olb — 2.6 
x 1(T 13 (10 4 K/T) 7 cm 3 s" 1 flKwan|l997| ). In this cal- 
culation, we have used olb — 3.9 x 1Q~| 3 cm 3 s _1 for the 
temperature of 5600 K (see |Kim fc Koo|200T ). N MV is the 
total number of ionizing photons per unit time, emitted 
by ionizing stars and "no" is the initial particle number 
density of the ambient neutral gas. We have adopted the 
Lyman continuum photon flux val ue (N w ; =) of 1.2 x 



10 49 ph s- 1 (logN n „ = 49.1) from |Kim fc Koo| fl2Q01| 
for an electron temperature, distance and integrated 21 
cm (1.43 GHz) flux density of 5600 K, 4.2 kpc, and 6.67 
Jy respectively. We have estimated the ambient density 
(n =) 3575.7 cm" 3 of the Hll region using 12 CO(J=3-2) 



line data (see subsection 3.2). Using N M? , , a mean radiu s 



of the Hll region (R =) 5^pc (see |Kim fc Koo||2001[ ) 
and n = 3575.7 cm -3 , we hav e obtained td yn ~ 3 . 3 Myr 
using Equation 1. Following, | Whitworth et al.| ([1994) 
analytical model for the "collect and collapse" process, 
we have estimated a fragmentation time scale (tf rag ) of 
0.86 - 1.73 Myr for a turbulent velocity (a s =) of 0.2 - 
0.6 km s _1 in the collected layer. We have found that 
dynamical age is larger than fragmentation time scale 
for no = 3575.7 cm -3 . We have plotted the variation of 
tfrag an d td yn with initial density (no) of the ambient 



neutral medium (see Fig. 12). The tf rag is also calcu- 
late d for different turbulent velocities in Whi tworth etl 
[ah] (fl994) model. It is to be noted that if td yn is larger 
than tf rag , then ambient density (no) should be larger 
than 900, 1500, 1700 and 1850 cm" 3 for different "a/ 
values of 0.2, 0.4, 0.5 and 0.6 km s _1 respectively (see 
Fig. [l2|) . We have also estimated the kinematical time 
scale of molecular bubbles of about 1 Myr (~ 4 pc/ 4 
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TABLE 2 

Physical parameters derived from SED modeling of the selected clustered YSOs (see the text). 



Source 


Age 






L* 


M env 


Mdisk 




Degeneracy/ 


Name 


log(yr) 


(M©) 


log(T(K)) 


log(L©) 


log(M yr- 1 ) 


log(M yr" 1 ) 


(mag) 


No. of models 


IRS 1 


5.17+0.72 


10.33+2.11 


4.23+0.26 


3.74+0.26 


-4.19+2.06 


-6.96+1.28 


19.55+5.47 


316 


IRS 2 


3.86+0.24 


22.40+2.61 


4.01+0.16 


4.49+0.20 


-2.96+0.12 


-5.97+0.51 


17.63+2.48 


158 


IRS 3 


4.92+0.37 


3.51+1.07 


3.66+0.03 


1.88+0.18 


-4.46+0.64 


-7.70+1.09 


13.19+2.08 


821 


IRS 4 


5.46+0.86 


6.75+2.53 


4.11+0.26 


3.13+0.56 


-4.93+1.84 


-8.77+2.71 


13.63+3.85 


40 


IRS 5 


4.50+1.09 


5.01+2.40 


3.79+0.27 


2.53+0.45 


-4.58+2.05 


-6.36+1.38 


19.49+5.13 


376 


IRS 6 


5.73+0.55 


2.47+0.97 


3.69+0.09 


1.24+0.29 


-6.15+1.62 


-8.49+1.68 


14.12+2.12 


3271 


IRS 7 


6.43+0.18 


4.51+0.75 


4.16+0.09 


2.52+0.24 


-7.90+0.17 


-12.63+0.83 


8.13+0.10 


154 


IRS 8 


4.60+0.59 


7.80+1.11 


3.80+0.17 


3.04+0.23 


-3.61+0.31 


-6.63+1.34 


18.00+5.82 


48 



km s _1 ), assuming bubble size of about 4 pc and veloc- 
ity dispers ion ~ 4 km s _1 from 12 CO(J=3-2) map (see 
subsection 3.2). The comparison of the dynamical age of 



the H II region and the kinematical time scale of expand- 
ing bubbles with the fragmentation time scale further 
supports triggered star formation and indicates the frag- 
mentation of the molecular materials into clumps due to 
"collect and collapse" process around the bubbles. Fur- 
ther, the spatial distribution and the clustering analy- 
sis of YSOs show the association of YSO clusters in the 
North, East, West and South-East regions with molecu- 
lar material along the bubbles. Among these, the aver- 
age age of sele cted 8 YSOs, including IRS 1 and IRS 2 
(see subsection 3.3.3), is about 0.12 Myr (estimated from 
SED modeling), which is less than the dynamical age 
of the Hll region. We have also found that these two 
sources (IRS 1 and IRS 2) are young and massive em- 
bedded protostars (about 10 and 22 M ) associated with 
the dense clump at the interface of the bubbles. It seems 
that the YSOs present on the edges of the bubbles are 
associated with the collected and fragmented molecular 
clum ps, possibly by the "collect and collapse" process 



clumps, possi bly by tne collect and colla r 
(see|Kang et al.||2QQ9| |Pomares etll]|2009l ' 
2011). This scenario is further supported b; 



Paron et al. 
►y the detec- 
tion of massive YSOs (IRS 1 and IRS 2) associated with 
the dense clump. However, one can not entirely rule out 
the possibility of triggered star formation by compres- 
sion of the pre-existing dense clumps by the shock wave. 
It is therefore difficult to distinguish the star formation 
scenario from the pre-existing condensations and/or by 
the "collect and collapse" process in our selected region 
with the avail able data pre s ented in this work. 

Recently, jOnaka et al.| (12009) reported very broad 
emission features around 22 /im in some MSF regions 
(viz., Cas A, Carina nebula and Sharpless 171) and sug- 
gested that this feature possibly originated from the dust 
grains formed in the supernova. Since the G8. 14+0.23 
region is about 0.3 degrees away from the W30 super- 
nova remnant (SNR) and Spitzer-IRS spectroscopic ob- 
servations are available for this region in the wavelength 
range of 20 - 36.5 /im, it is possible to explore the sig- 
nature of the supernova-induce d star formation act ivity 
around this region. Following |Onaka et al.| p009), we 
therefore re-examined the Spitzer-IRS processed archival 
spectra (AOR: 14928896; PI: Yoko Okada) obtained from 
"Spitzer Heritage Archive 3 " around the G8. 14 +0.23 re- 
gion (see IRS slit positions in Fig. [8] and also Oka da et] 
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Fig. 12. — The schematic diagram shows the variation of frag- 
mentation time scale (tf rag ) and dynamical time {tdyn) with initial 
density (no) of the ambient neutral medium. The fragmentation 
time scale is calculated for different turbulent velocity (a s =0.2, 
0.4, 0.5, &; 0.6 km s — 1 ). Dotted black lines indicate the value of 
"no" at which "t^n" is equal to u tf rag " for different "a s " values. 



see http:/ /irsa.ipac. caltech.edu/applications/Spitzer/SHA/ 



al. 2008). We however did not fin d any broad featur e 
around 22 /im (see also Table 2 of |Okada et al.| |2008). 
Hence, we may reject the possibility of the supernova m- 
duced star formation ac tivity in the G8. 14+0.23 region. 
Ojeda-May et al.| ( |2002D studied 9 MSF regions, includ- 
ing G8. 14+0.23, intne vicinity of the W30 SNR using 6 
cm radio observation and they also ruled out the possi- 
bility of a supernova induced star formation based on the 
6 cm detection around this region and age consideration 
of the W30 SNR. 

4. CONCLUSIONS 

We have explored the triggered star formation scenario 
around two bubbles associated with a southern massive 
star forming region G8. 14+0.23 using multi- wavelength 
observations. We find that there is a clear evidence of col- 
lected material (molecular and cold dust) along the bub- 
bles around the G8. 14+0.23 region. The JCMT CO ve- 
locity map reveals that the molecular gas in the bubbles 
is physically associated around the region G8. 14+0.23. 
The surface density of YSOs reveals ongoing star forma- 
tion and clustering of YSOs associated with the edge 
of the bubbles. We conclude that the YSOs are be- 
ing formed on the edge of the bubbles possibly by the 
expansion of the Hll region. We further investigated 
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the "collect-and-collapse" process for triggered star for- 
mation a round the G8. 14+0.23 r egion using analytical 



model of |Whitworth et al.l ( |1994[ ). We have found that 
the dynamical age (~ 3.3 Myr) ol the H II region is larger, 
and the kinematical time scale of bubbles (~ 1 Myr) is 
comparable to the fragmentation time scale (~ 0.86 - 1.73 
Myr) of accumulated gas layers in the region for 3575.7 
cm -3 ambient density. The comparison of the dynamical 
age and the kinematical time scale of expanding bubbles 
with the fragmentation time scale further supports trig- 
gered star formation and indicates the fragmentation of 
the molecular materials into clumps possibly due to the 
"collect and collapse" process around the G8. 14+0.23 re- 
gion. We have also investigated infrared counterparts of 
two young massive embedded protostars (about 10 and 
22 Mq) associated with dense clumps at the interface of 
the bubbles. It seems that the expansion of the Hll re- 
gion is also leading to the formation of these two young 
massive embedded YSOs in the G8. 14+0.23 region. Our 
study possibly favours the "collect and collapse" scenario 
for the formation of YSOs around the bubbles in the 



G8. 14+0.23 region. However, we can not entirely rule 
out the possibility of triggered star formation by com- 
pression of the pre-existing dense clumps by the shock 
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